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Langevin equation with a multiplicative stochastic force is considered. That force is uncorrelated, 
it has the Levy distribution and the power-law intensity. The Fokker-Planck equations, which 
correspond both to the Ito and Stratonovich interpretation, are presented. They are solved for 
the case without drift and for the harmonic oscillator potential. The variance is evaluated; it is 
always infinite for the Ito case whereas for the Stratonovich one it can be finite and rise with time 
slower that linearly, which indicates subdiffusion. Analytical results are compared with numerical 
simulations. 

PACS numbers: 02.50.Ey,05.40.Ca,05.40.Fb 

I. INTRODUCTION 

The Langevin formalism was introduced to describe motion of a particle which was subjected both to a Newtonian 
deterministic force and to the irregular influence of a bath of small molecules (the Brownian motion). That random 
force is uncorrelated and usually it is taken in the Gaussian form. Importance and wide applicability of the Gaussian 
distribution in the statistical physics follows from its stability: it constitutes an attractor in the functional space. 
According to the central limit theorem, a superposition of distributions with a finite variance leads to the Gaussian 
and its second moment is proportional to the time. However, many phenomena cannot be described in this way. 
Variance may depend on time stronger or weaker than linearly; as a consequence, the diffusion is anomalous 
Moreover, one can frequently encounter in nature systems far from thermal equilibrium, for which moments, in 
particular the variance, are divergent 0, S, i, H, Q ■ From that point of view, the Levy process, which is stable but 
may be characterised by divergent moments, is an important generalisation of the Gaussian process. 

Divergent moments of the stochastic driving force can be attributed to nonhomogeneous structure of the envi- 
ronment, e.g. fractal or multifractal, which produces long-range correlations. Levy statistics can emerge from the 
temporal nature of the underlying process due to a subordination to the ordinary Brownian motion, which is highly 
inhomogeneous The medium nonhomogeneity can enter the Langevin equation via the deterministic force. It is 
even possible to model random effects in this way (the quenched disorder P, Q)- However, there are also processes for 
which the fluctuations of the stochastic force directly depend on the state of the system. The autocatalytic chemical 
reaction, in which theproduction of a molecule is enhanced by the presence of the molecules of the same type, can 
serve as an example [9|. As a result, the fluctuating term in the Langevin equation is multiplied by some function 
of the macroscopic variables (the multiplicative noise). Some physical problems require taking into account both 
additive and multiplicative noise [llj . 

In this paper, we consider the Langevin equation with the multiplicative noise, 

x^F{x)+G{x)ri{t), (1) 

in which F{x) is the deterministic force and the uncorrelated stochastic force 'r]{t) possesses the Levy distribution. 
The case G(x) =const has been extensively studied. It has been demonstrated tl2| that both the force-free system 
and that driven by the linear force are described by the Levy distribution and then the variance is divergent. The 
strongly non-linear force, however, is able to confine Levy flights and make the variance finite Also problems 

with more complicated forms of F{x) were considered in context of Eq.lT]) with the Levy noise. They involve the 
barrier penetration and escape from the potential well [lH, as well as the transport in a Levy ratchet fl7|. The 
general case, with multiplicative noise, can be regarded as a result of the adiabatic elimination of fast variables for 
nonlinear processes with additive fluctuations. Eq.([T]) can also be directly applied to model specific phenomena with 
fluctuations which are characterised by long tails in the distribution and a variable, in particular power-law, intensity. 
It is the case, for example, in the field of finance, where Eq.((T]) could be a natural generalisation of the Black-Scholes 
equation 

In the general case of variable G{x) we encounter the well-known problem of interpretation of the noise term in 

Eq.([T]) and then of the Stieltjes integral J* G[x{T)]dri{T) which, in this form, is meaningless 07 HH- We must decide 
whether G[x(r)] is calculated before the noise acts, or after that. The former case corresponds to the Ito interpretation 

m 

/ G[x{T)]dv{T) - G[x{U^i)MU) - viU^i)], (2) 
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where the interval {0,t) has been divided in n subintervals (n —f oo). The above integral does not obey standard 
rules of the calculus, in particular the chain rule. Those rules are satisfied by the stochastic calculus introduced by 
Stratonovich. The stochastic integral [l^ in this interpretation includes the process value both at the beginning and 
at the end of each subinterval: 



/ G[x{T)]drj{r)^J2^ 



{T]{t,) - r]{U-{)\. (3) 



From the mathematical and technical point of view, the Ito interpretation is easier to apply e.g. in the perturbation 
theory [2^. It is a suitable choice if the noise consists of clearly separated pulses, e.g. for a continuous description of 
integer processes. If, on the other hand, a system has some finite correlations and the white noise is only an approxi- 
mation, the Stratonovich interpretation is more appropriate. It is the case if the noise has an external source, i.e. the 
noise source is not influenced by the system itself and it is possible, in principle, to turn off the noise [l^]. Possibility 
of applying standard rules of the calculus is an important technical advantage of the Stratonovich interpretation. It 
makes possible to solve the Langevin equation exactly for some nonlinear models with the multiplicative Gaussian 



white noise 25| . Physical predictions which follow from the Langevin equation with the noise in both interpretations 
can be qualitatively different. It is the case, for example, for the Ginzburg-Landau model with external multiplica- 
tive fluctuations, which describes noise-induced phase transitions caused by short-term instabilities of the disordered 
phase. The system exhibits that transition if one interprets the noise in the Stratonovich sense, but not if one uses 
the Ito interpretation ^2^ . We demonstrate in this paper that predictions of both formalisms are different also for 
the Levy diffusion process. 

The aim of this paper is to study the Langevin equation ([1]) with the multiplicative noise, which is given by the 
Levy distribution and the algebraic G{x). The Fokker-Planck equations (FPE) for both interpretations are presented 
in Sec. II. In Sec. Ill, the equation without external potential is solved and the case of the harmonic oscillator potential 
is discussed in Sec. IV. Results are compared with numerical simulations in Sec.V. 



II. FOKKER-PLANCK EQUATIONS 

In this section, we derive equations for the probability density distribution, which correspond to Eq.([T]) in both 
interpretations. The noise ?7(t) is the symmetric Levy stable distribution, defined by the following expression 

p^{x) = ^/2/^ exp(~K''k'')cos{kx)dk (4) 
Jo 

in terms of the order parameter /i (0 < /i < 2) and the generalised diffusion constant K > 0. If /i < 2, p,, has long 
tails: p,i ^ |a;|~''~^ for \x\ oo. The general Levy processes can be defined by the Lcvy-Khintchine formula which 
expresses the characteristic function in terms of the Levy measure v{x) [23] . In the symmetric and non-Gaussian case 
it reads 



Inp^(fc) 



{l-e'''=')iy{x)dx+ / (1 - e''''' + ikxy{x)dx 

\x\>l J\x\<l 



(5) 



and ^{x) — corresponds to the stable process. 

In the Ito interpretation, t he p robability density distribution is given by the fractional Fokker-Planck equation with 
variable diffusion coefficient [U 113 

§iPii^, t) = -^F{x)pj{x, t) + K^^^[\G{xTpi{x, t)l (6) 

where the Riesz-Weyl fractional operator [2^ is defined in terms of the Fourier transform: -^yr = 
Equations of the form ([6]) can describe also jumping processes. For example, the fractional equation with the variable 
diffusion coefficient follows from the master equation which models the thermal activation of particles within the 
folded, heterogeneous polymers [soj : variability of the diffusion coefficient results there from the intrinsic potential 
of the monomer. Moreover, Markovian versions of the continuous time random walk (CTRW) produce equations of 
the form ([5]) in the diffusion limit. In particular, the coupled CTRW model with a variable jumping rate, which 
describes Levy fiights in nonhomogeneous media, involves the fractional FPE ^Sl^l in the form . The diffusion term 
in Eq.® is, in this case, the jumping rate. The drift term may also appear if we allow for a non- vanishing mean of the 
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Levy distribution. The master equation describes also systems which are characterised by the internal noise. Those 
fluctuations emerge in systems of discrete particles and they are an inherent part of the very mechanism by which 
the state of the system evolves [24j . A precise form of the deterministic equation does not exist since it is impossible 
for systems with that noise to eliminate the fluctuations. Consequently, the master equation describes the evolution 
of the entire system stochastic process. 

The Stratonovich interpretation of the stochastic integral, Eq.Q, means that rules of the calculus ~ the chain 
rule and the ordinary variable transformation formula ~ can be applied. The stochastic variable can be determined 
by a stochastic equation with the additive noise which results from that with the multiplicative noise, obtained by a 
variable transformation. The above property of the stochastic integral can be proved on the assumption that the noise 
has a finite variance [l^, [2l|. If, in addition, trajectories are continuous (Lindeberg's condition), rules of the ordinary 
calculus apply to the Fisk- Stratonovich integral and one obtains a relation between integrals defined by Eq.® and 
([3]): they differ only by a simple additive term As regards the Levy stable processes, they can be approximated 
by processes with the finite variance by introducing a cut-off in the Levy measure iy(x)ui Eq.(l5]) (truncated Levy 
flights). Such an approximation is very accurate, also for a large value of the argument [32] ■ In Sec.V, we will present 
numerical examples which confirm applicability of variable transformation rules of the ordinary calculus for the Levy 
stable processes. 

Knowing the variable transformation rules, it is possible to transform Eq.([T]) to a new Langevin equation, which 
contains the additive noise, instead of the multiplicative one [ssj . For that purpose we introduce a new variable y 
and reduce Eq.((T]) to the form: 



where the transformation reads 



» = f (») + lit), (7) 







The corresponding FPE is of the form 



d d - 5'" 

^Ps(y, t) = -—F(y)ps{y. t) + Q^Psiv. t). (9) 

After solving the above equation, the solution of the original equation ^ follows from the probability conservation 
rule: 

ps{x,t)^ps{y{x),t)^. (10) 

For the Gaussian case, = 2, Eq.® can be easily expressed in terms of the original variable x and a direct 
relation between ps{x,t) and pi{x,t) can be established. The difference between Fokker-Planck equations for both 
interpretations resolves itself to the additional drift, K'^G{x)G' (x), called 'spurious' or 'noise-induced' drift. 

In the following, we solve Eq.(IT]) for two cases: without external potential and with the linear drift, on the assumption 
of both interpretations of the stochastic equation. We assume the diffusion coefficient in the algebraic form: 

G{x) = ixr''/^ (11) 

Results can be generalised to other forms of G{x) by applying the method from Ref. [s^l- 



III. FORCE-FREE CASE 



We begin with the case of the Ito interpretation. Eq.® with F{x) = becomes the fractional diffusion equation 
with the variable diffusion coefficient, 

dpi{x,t) ,,^ d'^[\x\-'pi{x,t)] 

-^^^^ ■ ('^^ 

The above equation results not only from the Langevin equation with the multiplicative noise. It constitutes the 
small wave number limil^ (the diffusion or fluid limit) of the master equation for a jumping process in the framework 
of the coupled CTRW 31]. That process is deflned in terms of two probability distributions: of the jumping size, in 
the Levy form, and the Poissonian, position-dependent waiting time distribution. The diffusion coefficient |G(a;)|^ 
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in Eq.Q is then the jumping rate and the parameter 9 governs the transport speed. In particular, for fj, = 2, when 
the variance exists, Eq. p^ describes the anomalous diffusion process: subdiffusion for 9 > 0, enhanced diffusion for 
9 < 0, and the normal diffusion for 9^0. The same classification holds also for /i < 2 when we introduce fractional 
moments [Slj]. 

Eq. can be solved in the diffusion limit by applying the Fox functions formalism [s^, [s^ [s^l • Details of the 
derivation are presented in Refs.jsil, [sH]. The Fox functions are well suited for problems which involve Levy processes 
since any Levy distribution, both symmetric and asymmetric, can be expressed as the function H2'2{x) [39j . Moreover, 
due to the multiplication rule, the term \x\~^pi{x, t) in Ea. (jl2p can be easily evaluated and it produces the Fox function 
of the same order, only the coefhcients are shifted. Those properties of Eq. (fT2|) suggest the method of solution: we 
assume the solution in the scaling form which involve the Fox function. 



p/(x,t) = Na{t)Hl 



a(t)|a 



{ai,Ai), (02,^2) 

(6l,i?l),(62,i?2) 



(13) 



where N is the normalization constant, and try to adjust the coefficients, as well as to derive the function a{t). To 
implement the approximation of small fc, we pass to the Fourier space, in which Eq. (|12p reads 



d ~ n ^ 



(14) 



According to the well-known formula, the Fourier transform from the Fox function is also the Fox function but of higher 
order: 2(0;)] = H^'l{k). Then we expand the Fox functions, which correspond to pi{k,t) and Tc[\x\~^pi{x,t)], 

in the fractional powers of k; terms of the order \k\'^^~^^ and higher are neglected. Both sides of Ea. p4| can be 
adjusted only if all terms except k'^ and |fc|'^ vanish. We can eliminate adverse terms by a proper choice of the Fox 
function coefhcients. Inserting the coefficient values, determined in this way, to Eq.(fT2|). yields the solution in the 
form 



pi{x,t) 



Ma{t)Hlj 



a{t)\x\ 



(1 



t2+0 ■ tj.+e 
(61, Si), (1 - 



(15) 



the initial condition is pi{x,0) = 6{x). Moreover we assume fi + 9 > 0. The coefficients (02,^2) and {bi,Bi) are 
essentially arbitrary and one needs additional requirements to settle them. Expansion of the Fox function in Eq. (flT 



yields pi{x,t) ^ \x\''^^^^ {x 0). It explains why the coefficients (6i,i?i) have not been determined: since small 
fc| corresponds to large the diffusion approximation does not cover the region of small In the following, we 
assume bi = 9 and Bi — 1, which choice is exact for CTRW [11]. The coefficients (02,^2), in turn, correspond to the 
shape of the distribution in the limit /i — > 2 [s^. From Eg. p^ follows a simple differential equation for a{t) which 
yields a{t) ^ ^^1/(^+6), xhe asymptotic expansion of Eq. (fT5|) yields pi{x,t) f^^^f^^^^x]^^^^ . Therefore, we obtain 
the Levy process which has the same distribution as the driving process, Eq.Q. The variance and all higher moments 
are divergent. 

For the Stratonovich interpretation, we introduce the new variable y, according to Eq.®, 



l2.|(f+e)/Msgna;, 



if (/i + 9) ' 

which ends in Langevin equation with the additive noise. Then Eq.® reads 



It can be solved exactly and the solution expressed in the form [40|, |41| 



psiy,t) = -^H^j 

M|y| ' 



(1,1^), (1,1/2) 
(1,1), (1,1/2) 



(16) 



(17) 



(18) 



which corresponds to the symmetric Levy process y{t). The inverse transformation to the original variable yields 



^^\xn-' 



{l + 9/fi){KHy/t' 



(1,1/m), (1,1/2) 
(1,1), (1,1/2) 



(19) 
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The asymptotic expansion of the Fox function in Eq. fTO)) is given by the expression ps{x,t) ~ t'^/(t^+^')\x\~^~'^~^ 
{\x\ oo). It differs considerably from the Ito result: the shape of the tail depends not only on the order parameter 
of the driving process, /x, but also the 0-dependence emerges. As a result, the variance may not be divergent. We can 
express the variance by Mellin transform of the Fox function, x{-^)i ^J^d evaluate it by simple algebra: 



(x^) =2 / x^p{x,t)dx = 2 
Jo 

- k(- + 1) 



K{- + 1] 



2fi 



(20) 



2 

TT/J, 



Ki- + 1) 



2m/(m+0) 



r 1 



2/i 



where we assumed 5 — Iji/lji + 9) < fj,, which implies 6 > 2 — fj,. On that condition, the integral in Eq. pD| is 
convergent and the variance exists. We conclude that the diffusion process - in which the stochastic driving force is 
Levy distributed and the stochastic equation is interpreted in the Stratonovich sense - may not be accelerated for the 
case without any external potential. If the variance exists, the diffusion is anomalously weak since the convergence 
condition coincides with the subdiffusion condition: the variance rises with time slower than linearly, {x"^) ^ . The 



slope of that dependence diminishes with the parameter 9. In the case /i 
enhanced diffusion occurs, for 6* < 0, as well as the normal one if ^? = 0. 



2, beside the subdiffusion, also the 



IV. LINEAR FORCE 



In the case of stochastic motion in the harmonic oscillator field, F(x) ~ —Xx, which is governed by the Langevin 
equation with the Gaussian white noise (the Ornstein-Uhlenbeck process), the probability distribution converges with 
time to a steady state which corresponds to the Boltzmann equilibrium distribution. If the driving noise has the Levy 
distribution with /i < 2, the stationary limiting distribution still exist but the Boltzmann equilibrium is not reached 
and the variance is infinite [T2| . We will demonstrate that, if the multiplicative noise is interpreted in the Stratonovich 
sense, the steady state may have the finite variance. 

In the Ito interpretation, FPE is given by Eq.®, 

r) rl Bt^ 

-pj{x, t) = A-[xpK2^, t)] + V(:r, t)l (21) 

and its Fourier transformation yields 

-Piik, t) = -Xk—piik, t) - K>^\k\'^T,[\x\-%i{x, t)]. (22) 

The solution of Eq. ip^ can be obtained ^4§| in a similar way as in the case F{x) = 0, namely by inserting the 
expression (I13p into Ea. (j22p . Then the Fourier transforms of the respective functions are expanded and terms of the 
order \k\^^^^ and higher are neglected. The same conditions for the Fox function coefficients are required because 
contribution from the drift term contains only the component Therefore we obtain the solution of Ea. (PT|) in the 
form p^ . The comparison of terms of the order \k\^ on both sides of Eq. (|^^ results in a simple differential equation 
for the function a{t). Its solution reads 



a{t) 



X/CL 



1 - cxp[-A(/x + 9)t] 



(23) 



The constant cl — K^^ho/ iih^^ involves the expansion coefficients /i^ and of the functions pi and .Fc[|x| ^pi], which 
correspond to the orders \k\^ and respectively. They are given by /i^ = N{iJ,+9)T{—ij)r{l+iJ.+9) cos(7r/x/2)/r[a2 + 
^2(1 + fi)]r[-{^i + 9)/{2 + 9)] and ho = N{fi + 9)/{2 + 9)T[a2 + ^2(1 - 9% where the normalization constant 
N = r[-6'/(2 + e')]r[a2 + A2]/2r(l + 9)T[-9/{ii + 9)]. The numerical values of the solution can be computed by 
means of series expansions, both for small and large Expansion of the function (jlSp in powers of \x\~^ produces 
the following expression 

't{n^2 + A2[\-9+{^i + 9)^])T{->^^)^\^' 



6 



10 



X 

V 4 







1 1 1 1 1 1 1 
1-1=0.5, X-2 






H=0.5, X-^ 






1.1=0.5, ?.=0.1 






H=1.5, /.=2 






- - n=1 .5, X=-\ 


- I 




H=1.5, X=0.1 








1 . 1 . 1 . 1 . 



10 



FIG. 1: Variance of the steady-state solution of Eq.([T| in the Stratonovich interpretation for the case F{x) 
from Eq. pO)) . 



'Xx, calculated 



'a. 




FIG. 2: (Colour online) Comparison of numerical solutions of Eq.(IT} for the case G{x) — x and F{x) = 0, with fi = 1.5 at 
t — 1. The distribution pi{x,t) (lower line) was calculated according to Eg. (131) and ps{x,t) from Ea. p2)l (upper line). The 
dashed line presents result of the evaluation of ps{x,t) by means of Eq.(|3j. 

where a{t) is given by Eq. ()23|) . In the following, we assume the remaining Fox function coefficients in the form: 
02 — 1/2 + 9{\ + 9)/{2 + 6) and A2 = 1 — 1/(2 + 6). For these coefficients, pi{x,t) agrees with the exact diffusion 
equation sohition in the case fi = 2: it corresponds to the stretched-Gaussian dependence [ssj 



pi{x, t) - a^+^\x\^ exp[-const(a|a;|)2+'^] 



(25) 



Eq. ([24|l implies the asymptotic shape of the distribution, pi{x,t) ~ |a;|~^~^, the same as for the case F{x) — 0. 
Therefore, the variance is also divergent for all 6 and ^ < 2. 

To obtain the solution of Eq.([T]) in the Stratonovich interpretation, we introduce the variable y, which is given by 
Eq. (fT6|) ■ and transform the drift term according to Eq.®. The Langevin equation takes the form 



A 



(26) 
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FIG. 3: (Colour online) Comparison of numerical solutions of Eq.(IT} for G{x) — Ix]'^^'^ and F(x) = —Xx, with parameters: 
fi = 1.8, A = 1 and t = 1, in the Stratonovich interpretation. Distributions marked by lines were calculated by using the variable 
transformation, from Eq. (|32|) . and those marked by dots - directly from Eq.Q. Cases for two values of 9 are presented: 9 — —0.2 
(upper line for small and large a;) and 6 = 0.5. 



and FPE, expressed by the new variable, has the constant diffusion coefficient: 

-ps(y,i) = A(l + e/f,) — [ypsiy.t)] + ^ps{y,t). 



(27) 



The above equation can be solved exactly with the initial condition psiy^ 0) — 6(jj) The Fourier transform of the 
solution reads ps(k,t) — exp(— i(r^(T(t)|fc|^), where 



1 - exp[-A(/u + e)t] 
A(/i + 9) ■ 



(28) 



After inverting the Fourier transform and transforming back to the original variable, we obtain the probability density 
distribution in the following form 



1,1 

21 I "-"-2.2 

p. \x\ 



-H. 



ii+e/M 



(1 + 0/fi)Ka^/t^ 



(l,l//i), (1,1/2) 
(1,1), (1,1/2) 



(29) 



which can be evaluated by series expansions, similar to Eq. pi]) . The tail of the distribution ps{x,t) has the same 
form as in the case F{x) = 0: ps{x,t) |a;|^^^^~^. The second moment is convergent ii 9 > 2 — /i. On that condition, 
the system reaches with time a steady state which is characterised by the variance 



lim {x'^){t) = 



K{- + 1) 



TTfJ, I p 



2^i/(^i+e) 



r 1 



2p 



np 

p + t 



[i^l + 9)X]-^^<^'^+'l 



(30) 



The above quantity is presented in Fig.l as a function of 9 for some values of p and A. In all cases it declines with 6; 
this fall is particularly rapid for large A. Predominantly, the distribution is broader for smaller p but this trend turns 
to the opposite in the limit ^ oo. The parameter 6 influences the convergence speed to the steady state, according 
to Eg. ([^5]) . which is the same as for the Ito interpretation, cf. Eq. 



NUMERICAL SIMULATIONS 



In this section, we compare the analytical results with numerical stochastic trajectory simulations from the Langevin 
equation ([T]). For the Ito interpretation we apply the Euler method. Eq.([2]) implies that for each integration step 
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1 10 
X 



FIG. 4: (Colour online) Probability distributions obtained from trajectory simulations for the case F{x) = — Aa;. The 
Stratonovich {ps{x,t)) and Ito (pi{x,i)) results are marked by the dashed and dashed-dot lines, respectively. The analyt- 
ical solutions of Eq.lTJ (solid lines) follow from the series expansions of Eq. (|29|) both for small and large |a;|, as well as from 
Ea. (|24|l . Parameters are the following: ^ = 1.2, 6 = 1, X = 1, and £ = 1. 

T the stochastic integral from the function G(x) can be expressed by the noise value rji by means of the following 
formula [H 

l-t i-t N 

/ G{x{t))f^dt^ / G(a:(^))d7, = ^G(x(^,))Tl/''^7^, (31) 
Jo Jo 

where ti = it and N = t/r; the random numbers rji are sampled from the Levy distribution ^ with the order 
parameter /i. In the case of the Stratonovich interpretation, we first transform Eq.([T]) to the corresponding equation 
with the additive noise. Then we apply the Hcun method of integration, 

= + [F{y{U)) + F[v{U) + F{y{U))T]]T /2 + t^'^^ij,. (32) 

Transformation back to the original variable produces the stochastic trajectory x{t). We demonstrate results of those 
algorithms in Fig. 2 for the simple case G{x) — x and F{x) = 0. The applied parameters, number of trajectories 10^ 
and T — 0.001, ensure sufficient accuracy. The above results are compared with the probability distribution obtained 
from the numerical integration of the equation with the multiplicative noise, in which the stochastic integral is defined 
by Eq.(I3]). Agreement of both results is very good. A similar comparison is presented in Fig. 3 for the nonlinear G{x) 
in the form (|lip with two values of the parameter 9, both positive and negative; no cut-off in the Levy measure was 
introduced. This case is numerically more complicated because the difference formula is not explicit and numerical 
solving of a nonlinear equation is required at each integration step. Agreement of both methods of calculation in 
Fig. 3 demonstrates that rules of the ordinary calculus are applicable for the Levy processes with infinite variance. 

The comparison with analytical results for the linear drift, presented in Fig. 4, indicate a good agreement of those 
methods of solution. On the other hand, results which correspond to the Ito and Stratonovich interpretations for 
the same case differ considerably. Since they coincide for 9 = 0, one can expect that the difference rises with \9\. 
In Fig. 5, the probability distributions for various 9, evaluated by numerical trajectory simulations, are presented. 
The difference between pi{x,t) and ps{x,t) for 9 = 2 is indeed very large. The slope of pi{x,t) remains constant 
for a given /i and that of ps{x,t) rises with 9. We present those slopes, as a function of the order parameter /i, in 
Fig. 6. The slopes rise with /i, according to the analytical results — /i — 1 and — /i — 6* — 1 for the Ito and Stratonovich 
interpretations, respectively. 

VI. SUMMARY AND CONCLUSIONS 

We have studied the non-linear Langevin equation with the multiplicative white noise which is distributed according 
to the Levy statistics and has the power-law strength. The corresponding Fokker-Planck equation is fractional and its 
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FIG. 5; (Colour online) Probability distributions ps(x,t) and pi{x,t) obtained from trajectory simulations for the case F{x) = 
—Ax at t = 1, with /i — 1.5 and A = 1. 




FIG. 6: Slope of the solution pi{x,t) (lower points) and of ps{x,t) (upper points) for large |a::|, calculated by the numerical 
integration of Eq.([l} for the case F{x) — ~Xx and t = 1. The parameters: A = 1 and 9—1. The lower and upper straight 
lines have the slopes — /i — 1 and — /i — 6 — 1, respectively. 



form depends on the interpretation of the stochastic equation. In the Ito interpretation, FPE possesses the variable 
diffusion coefficient. In the case without any external force, FPE coincides with the diffusion equation which was 
obtained in the framework of the coupled CTRW with the position-dependent waiting time. Solution of FPE with 
variable diffusion coefhcient, both for the case without drift and for the harmonic oscillator potential, represents the 
Levy process with simple scaling and the same order parameter as the driving noise. This property does not hold 
for other potentials. For example, solution to the problem of the wedge-shaped potential, F(x) sgnx, studied in 
Ref. [131, is a combination of two Levy processes and simple scaling is lacking. Nevertheless, slowly decaying Levy 
tails are present in that solution and then the variance is divergent for all 6. 

In the Stratonovich interpretation, FPE has been obtained by the variable transformation; in this case the problem 
is reduced to Langevin equation with the additive noise. The resulting probability distributions differ considerably 
from those in the Ito sense. In particular, shape of the tail depends on the parameter of noise intensity 6 and, as 
a result, the variance may be convergent. Therefore the diffusion process, for the case without drift, can be either 
accelerated or anomalously weak, in contrast to the Ito result, which predicts only the accelerated diffusion. For the 



10 



case ^ = 2, both approaches differ by an additional, effective force in the Langevin equation. The disagreement between 
both interpretations is in fact not surprising since the deterministic force in the Langevin equation is not just the 
Newtonian one; those forces are identical only in the case of the external noise, i.e. for the Stratonovich interpretation 
[2?|. Since the difference between the interpretations is deterministic in nature no different underlying stochastic 
dynamics is implied. For the general Levy stable processes, relation between distributions in both interpretations is 
more complicated and the corresponding equations cannot be related one to the other by means of a drift term. 

The distribution tails, which are algebraic and fall rapidly enough to ensure convergence of the variance ('fat tails'), 
are of physical importance. They are well known e.g. in the field of the economic research ^44 . ,45] . Based on the 
poor empirical performance of the Black-Scholes model of option pricing, which mathematically resolves itself to the 
Langevin equation, one proposes to replace the Gaussian noise by the a Levy one but with a truncated tail [lB|- Such 
process may converge to the Gaussian so slowly, that numerical calculations yield only fat tails [H, . The Langevin 
equation with the multiplicative Levy noise in the Stratonovich interpretation could be an alternative model of the 
fat tails. 
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